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Abstract 

We derive the spectral dependence of the non-linear susceptibility of any order, generalizing 
the common form of Sellmeier equations. This dependence is fully defined by the knowledge of 
the linear dispersion of the medium. This finding generalizes the Miller formula to any order of 
non-linearity. In the frequency-degenerate case, it yields the spectral dependence of non-linear 
refractive indices of arbitrary order. 
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I. INTRODUCTION 



Non-linear optics [l| relies on the knowledge of the non-linear susceptibilities (or, alter- 
natively, the non-linear indices) of the propagation media. This description is generally 
truncated to the first term, i.e. the second-order susceptibility in non-centrosymmetric me- 
dia, or the third-order susceptibility in centrosymmetric ones. However, the increase of the 
available laser powers as well as the investigation of systems like optical fibers [2] or pho- 
tonic crystals [3| where the confinement of the light increases its intensity raise the need to 
consider higher-order processes. Recently, the non-linear refractive indices of O2, N 2 were 
measured up to n 8 (i.e. the 9 i/!, -order susceptibility x*- 9 '*) an d in Ar up to n w (i.e. x 



jp 



(11) > 



rat they must be considered in the de- 



10] . This result, which was unexpected, 



at 800 nm |4j. Furthermore, we demonstrated |5| 
scription of the filamentation of ultrashort pulses 6- 
provides a clear illustration of this need, and of the associated requirement to evaluate these 
terms at any wavelength. But systematic measurements over the spectrum are out of reach 
of present experimental capabilities. A theoretical support is therefore required to extend 
the existing experimental results to any frequency, which would provide a new insight into 
non-linear optics. 

Such relation was provided by Miller in the case of the second-order susceptibility. He 
observed that in many crystals the knowledge of x f° r one single triplet of frequencies 
(tu' ; u[, co' 2 ) and the dispersion relation of the medium is sufficient to determine x^ f° r an y 
triplet of frequencies (tuo; cox, CO2) [llj : 

x (2) (^o;^i,^2) _ x (1) (^o)x (1) (^i)x (1) (^2) m 
x^H-u' x ,u' 2 ) - x (1) K)x (1) K)x (1) K) 1 } 

This formula has been widely used e.g. determine the to derive the second order static 
susceptibility in semiconductors 12(. Mizrahi and Shelton later suggested that the Miller 
formula could be extended to provide the spectral dependence of the non-linear index n,2, i.e. 



to the third-order frequency-degenerate non-linearity [13|. However, they did not demon- 
strate this result in their article, but rather referred to an unpublished work [14]. Later, 
Bassani et al. demonstrated this relation for any order of non-linearity in the specific case 



of harmonic generation processes [15 
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In this Letter, provide an explicit expression of the spectral dependence of the non-linear 
susceptibility of any order, generalizing the common form of Sellmeier equations. We show 
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that this spectral dependence is fully defined by the knowledge of the linear dispersion of 
the medium. As a consequence, the Miller formula ([I]) ll| can be extended to any order of 
non-linearity. In the frequency-degenerate case, this finding yields the spectral dependence 
of ion-linear refractive indices of arbitrary order, confirming and widely generalizing Mizrahi 
and Shelton's statement about x*- 3 '* Lli 



II. DERIVATION OF GENERALIZED MILLER FORMULAE 



Electrons in an anharmonic oscillator 



Elucidating the spectral dependence of the electric susceptibility requires to consider the 
equation of motion of an electron located at Fin a three-dimensional potential V(r). V can 
be expanded as a 3-dimensional Taylor series around the equilibrium position r = 0: 



v(f) = v(6)+ Yl 



X lyJ Z k 

i\j\k\ 



dx l dyidz k 



(2) 



J r=0 



dV 
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y=0 



[f] 
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i+j+k>2 

where the summation begins &tq = i + j + k = 2 since x _ Q 
by definition of the equilibrium position. As a consequence, the electron experiences a force 
equal to 



F(r) = -W(r) = -J2 

q>2 



( x 1 1 y J z k ^ 
(i-iy.jlkl 



! 1p 1 Z k 



d q V 



dx l dyidz k 



f=0 



i\(j-iy.k\ 
\i\j\(k-iy./ 

The macroscopic polarization along the axes x, y and z are respectively P x 



(3) 



—Nex, 



P y = —Ney and P z = —Nez, N being the local density of electrons and — e their charge. 



Equation ([3]) therefore rewrites: 



q>2 



or, equivalently: 



Ne) 
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d q V 
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dx l dyi +1 dz k 

di +1 V 
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(5) 
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We introduce the (q + l) th rank tensor Q^ q \ which elements are given by: 



^ x;x<j) , y U) , Z W 

{q) 



/ 



m x z!j!fc! 

\ dx l dyi dz k + 1 jt—fi/ 

where, for example, indicates that the coordinate x appears i times in the index list of 



di +1 V 



dx l+1 dyi dz k 

di +1 V 
dx l dyi +1 dz k 

di +1 V 



\ 



J r=0 
r=0 
r=0 



(6) 



the considered tensor element. As a consequence, the classical equation of motion of the 
electron becomes: 



H 2 P HP 00 q P 



dt 2 



9=1 



1 = 1 



iVe 2 



m 



E(t) 



(7) 



field, u 2 



<v 



\ 



is the eigenfrequency matrix of the considered medium, which 



where : denotes the contracted product and (^) the tensorial product. E is the driving electric 
/ , ..2 n ^ 



ul z j 

is diagonal provided y, and 2; are the principal axes of the optical frame of the medium 
(e.g. u e>x = J ^fj^r), and the matrix 7 stands for the linear absorption. No multiphoton 
absorption is considered here. Note that a medium symmetry of at least KXC2 (te. a 
crystal with at least orthorhombic symmetry, or a statistically isotropic medium like air) 
allows to simultaneously diagonalize Co 2 and 7 in the optical frame 17]. 
A perturbative solution of equation (J2j) is searched as: 



J2 <x l P {l) 



(8) 



1=1 



where a e]0..1[ is a free parameter. The series begins at 1, assuming = 0. Inserting it 
in equation (JTj) and equating the terms in a q yields a set of equations 



d 2p(i) _ d p(i) iVe 2 -. 

+ 7—^ + w e 2 P (1) = —E(t) 



dt 2 



dt 



m 



and, Vg > 2, 



^— + 7^— + cD e 2 P (9) = -NeQ^ 



dt 2 



dt 



u=l 



p(l) 

~Ne 



(9) 



(10) 



Note that in the right-hand side of Equation (fTUj) . we have deliberately omitted the 
terms in ®j] ? '=q (-P^), with g > > 1. These terms imply non-linear polarizations 
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in the construction of the considered higher-order non-linear polarization and hence cor- 
respond to cascades of frequency mixings, like e.g. in the generation of third harmonic 
by frequency-doubling the fundamental wavelength and then mixing it with its second- 
harmonic. Consistently, in the identification of the non-linear susceptibilities in Equation 
(12Uj) . only single-step mixing will be taken into account. 



B. Sellmeier equations 

If Cl can be diagonalized (i.e. if the medium symmetry is at least KXC2 or absorption can 
be neglected), Equation (Q defining the linear polarization is purely vectorial and can be 
solved on each axis independently. Omitting the indices displaying the axis of the considered 
polarization, we write PW(i) = J^°° P^ l \uj)e %u)t duj and E(t) = J^°° Eo^e^doo, where, fol- 
lowing the usual notation, positive frequencies denote incident ones and negative frequencies 
are emitted ones. Introducing the spectral dependency parameter Q(u) = uj 2 — u 2 + iuj 
yields 

r+00 at- 2 r+oa 

/ n(u)P^\io)e luJt duj = / E (u)e Mt du (11) 

J— 00 m J -00 

which implies, for any frequency Uq of the emitted field except for the resonance frequency : 
Hence, the linear susceptibility at frequency coo is: 

x () M = Try — r — = — r~ 2 — t—- — V 13 

which defines the spectral dependence of x 011 an y principal axis. If absorption can be 
neglected (7 ~ 0) and x ^ 1; Equation (fl3|) takes a form similar to a typical Sellmeier 
formula: n 2 — 1 = gzfrp ; where n 2 = 1 + x' '• 



C. Uniaxial generalized Miller formula 

Equation (flOj) can be reduced to a scalar form provided the polarization is excited along 
one single principal axis. We will first focus on this case which simplifies the writing and 
thus helps focusing the discussion on the the physical aspects, without altering the principle 
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of the derivation. Omitting the index corresponding to the axis, the scalar form of Equation 
ffTUj) reads, for any q > 2: 

1 y + 7 ^±_ + ^ P ii) = -NeQ® l^—) (14) 



dt 2 ' dt 



Since P^(t) = /_ ~ P w (w)e^*rfw, Equation (ED rewrites: 



(iVe) 



i //-//n^Me^) (16) 



/=i 

Identifying the terms at an arbitrary frequency Wo on both sides of the equation, we 
obtain: 

(3^E„(n«) 

Let us now consider the construction of a wave at Co>o from a set of q' < q incident waves 
at frequencies ui,...,u q i, each algebric frequency ui being implied times in the process: 
YH=i u i = Q> while energy conservation imposes YH=i u i u i = u o- In this case, 

-*n(Hfy 

where Cq 1 ' "'"' = q\/{u\\ x .. x u q i\) is the number of combinations achievable from q' sets of 
Mi, ..,ty objects, respectively. The terms of Equation ( TT~8]) corresponding to each combina- 
tion of frequencies conserving energy can be identified with the expression of the non-linear 
polarization using the (q + l) th rank tensor. 

P (9 Vo) = eoCT'--' % 'x {q) ^o;co 1 ,...,co q )l[E^(co l ) (20) 

i=i 

As stated above, consistently with Equation ( fl4|) we do not consider here the cascaded 
processes. Therefore, the identification yields: 

x"W-.«.) = - M ^nf'L,) pi) 



me 



o 



Niei +1 



Q (9) n^ (i) (^) ( 22 ) 



This expression provides a general description of the non-linear susceptibility of any order, 
provided the shape of the potential V is known. In practice, this potential is rarely known, 
but it is independent from the excitation frequency. The spectral dependence of the above 
expression is therefore only driven by Q, i.e. the spectral dependence of the first-order 
susceptibilities. As a consequence, in a given medium, the knowledge of both the frequency 
dependence of the linear susceptibility and of the q th -order susceptibility for a specific set 
of wavelengths is sufficient to extrapolate this susceptibility to any other set of wavelengths, 
through the relation: 



When applied to y_ , this generalized Miller formula immediately reduces to the original 
one of Equation ([I]) [11]. 

D. Three-dimensionnal Miller formulae 

While the one-dimensional treatment of Equations (fl"4"]) -f l23"]) provides an easy writing 
of the derivation, the same sequence can be applied to solve Equation fTTOT) in the general 
three-dimensional case. Products simply have to be transposed to the adequate tensorial 
products. The identification of the spectral components of the Fourier expression of 
and (£) yields a three-dimensional equivalent of Equation ( TT8l) : 



that this expression can easily be identified with the counterpart of Equation (1201) expressing 
the non-linear polarization in terms of non-linear susceptibility: 



(23) 




(24) 



In media with at least KxC2 symmetry, or if absorption can be neglected, f2 is diagonal so 




:(g)£(^) 



(25) 



i=i 



n?=o 



(26) 



(27) 
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where each index vi denotes either the x, y, or z axis. The spectral dependence of x 
only depends on the product of the Q Vl {ui), meaning that all frequencies commute. This 
property, together with the fact that the fl{co) are even functions as soon as absorption can 
be neglected, provides a direct evidence and generalization to arbitrary orders of the ABDP 
relations. These relations state that as soon as Cl can be diagonalized, Uq commutes with 



any uoi in the expression of x^ an d X^ 18 j . 



E. Application to non-linear refractive indices 

In the case of frequency-degenerate interactions of odd order excited along one single prin- 
cipal axis, the susceptibility tensor reduces to x^^K^o) = ^ (2p+i)| 1 ' >! °no°(io) n 2 P (^o); 
which defines the p* h -order non-linear refractive index ni p . Therefore, the knowledge of the 
dispersion curve for the linear susceptiblity and the measurement of x^ 2p+1 ^ at one single 
frequency u' are sufficient to provide x at any frequency ojq thanks to Equation (1231) : 

or equivalently, if absorption can be neglected, in terms of the real parts of the non-linear 
refractive indices: 

n 2p {u ) _ ( ^o(^o) - l\ 2p+2 ^ 



n 2p (u' ) V^oK) - 1 

As an illustration, Figure [1] displays the values of n 2 through n% for N2,Oo and Ar, based 
on the recent experimental measurements of Loriot et al. at 800 nm |4|, |5j, extrapolated 
to the whole spectrum using Equation (1291) and the dispersion data of Zhang et al. [3]. 
These dispersion data have been chosen because they have been measured consistently for 
all proposed species. However, we checked that they agree within 2% with typical data 
available over the spectral range plotted Q-Q. 



III. GENERALIZATION TO MIXES AND MULTIPLE RESONANCE FREQUEN- 
CIES 

Up to now, we have only considered the case of a single type of oscillators. However, 
in actual media {e.g. in mixes like the air, where both N 2 and 2 get polarized, or in 
crystals where several oscillation modes can be excited) several eigenfrequencies uj e ^ may 
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FIG. 1: Spectral dependence of the non-linear refractive indices (a) ri2, (b) n 4 , (c) n 6 , and (d) rig 
of O2, N2, air and Ar at atmospheric pressure 

— * — * 

contribute. In this case, the polarization P can be defined as 22 Pi, where I denotes the 
types of oscillators, or species. If the coupling between the oscillators can be neglected, the 
above derivation applies to each oscillator type independently, and the resulting refractive 
index will be given by the Lorentz-Lorenz model. For example, away from the resonances, 
the refractive index of air is given by the Sellmeier equation |l9j |: 

in8 , ^ onirr1 , 2368616 19085.73 

10 8 (n - 1) = 8015.514+ — + — (30) 

vo ' 128.7459 - 1/A 2 50.01974 - 1/A 2 V ; 

where the wavelength A is expressed in fim. Obviously, the terms of Equation ( 1301) respec- 
tively correspond to N 2 , with a resonance at 1/A 2 ~ 128.7/im -2 (i.e. A = 88 nm): 

H/ s 2398095.2 . . 

10 8 (n , N2 - 1) = 8736.28 + (31) 

v2 _ , cn,,™-2 



and to 2 , with 1/A 2 ~ 50/im- 2 (i.e. A = 141 nm) [19[: 

H/ s 456402.97 . , 

10 8 K, O2 - 1) = 15532.45 + — _ 2 (32) 

In such cases, the generalized Miller formula of Equation f l23|) cannot be applied to the 
material or the mix as a whole. Instead, they must be applied to each susceptibility-order 
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of each oscillator type. The non-linear susceptibilities of the whole material will then be 
deduced from those of the individual oscillator types through the Lorentz-Lorenz model. 
Figure [TJ displays the result of this treatment in the case of air. 



IV. CONCLUSION 

In conclusion, we have explicitly derived a generalization of the common form of Sellmeier 
equations, in both frequency-degenerate and non-degenerate systems. This generalization 
provides the spectral dependence of the non-linear susceptibility of any order, which is 
fully defined by the knowledge of the linear dispersion of the medium. As a consequence, 
the Miller formula ([1]) can be generalized to any order of non-linearity and any tensor 
element of non-linear susceptibilities as soon as the material has at least KXC2 symmetry and 
negligible absorption. In particular, the spectral dependence of non-linear refractive indices 
of any order can be obtained from their value at one single frequency and the dispersion 

lying 



curve of the medium. Such knowledge is of particular value in non 
confined light leadin 
or photonic crystals 



mear optics imp 

confined light leading to very-high intensities, as, e.g. in fiber optics [2] fllamentation [7 
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